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Abstract 

It has been experimentally established in numerous cases that precip- 
itation of mono dispersed colloids from homogeneous solutions is a com- 
plex process. Specifically, it was found that in many systems nuclei, 
produced rapidly in a supersaturated solution, grow to nanosize primary 
particles (singlets), which then coagulate to form much larger final col- 
loids in a process dominated by irreversible capture of these singlets. This 
paper describes a kinetic model that explains the formation of disper- 
sions of narrow size distribution in such systems. Numerical simulations 
of the kinetic equations, with experimental model parameter values, are 
reported. The model was tested for a system involving formation of uni- 
form spherical gold particles by reduction of auric chloride in aqueous 
solutions. The calculated average size, the width of the particle size dis- 
tribution, and the time scale of the process, agreed reasonably well with 
the experimental values. 
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1. Introduction 

Ever since Faraday reported in 1857 the preparation of gold sols 
of different colors (1), scientists have been fascinated by monodispersed 
colloids. However, systematic studies of the synthesis, properties, and 
mechanisms of formation of such colloids have been initiated only quarter 
a century ago. Since then, a large number of uniform dispersions of 
particles of simple and mixed chemical compositions and various shapes, 
ranging in modal size from several nanometers to several micrometers, 
have been described in the literature. The method of choice has been 
precipitation from homogeneous solutions, either directly or via gel/sol 
or sol/gel routes (2-7). 

Two essential ingredients are necessary to explain formation of uni- 
form colloids by direct precipitation. The first is to establish the identity 
of the solute species that are involved in the solid phase formation. The 
second is to explain all stages leading from the initial nucleation to fi- 
nal particles. Early works on the precipitation of monodispersed colloids 
were based on the concept suggested by LaMer, which essentially assumes 
a short nucleation burst, followed by diffusional growth of the nuclei to 
form identical, finely dispersed particles (8,9). This mechanism appeared 
simple and plausible and was generally accepted as operational. 

With the increasing availability of monodispersed colloids of spheri- 
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cal and other shapes, it has become evident that the burst-nucleation/dif- 
fusional growth mechanism alone cannot explain many experimental ob- 
servations. The first among these has been the finding that some per- 
fectly spherical particles, directly precipitated from solution, showed X- 
ray characteristics of known minerals, such as ZnS of sphalerite (10), CdS 
of greenockite (11), Fe203 of hematite (12), etc. It was obvious that these 
particles could not be single crystals and, indeed, it has been confirmed 
by different techniques that a great majority of monodispersed colloids 
consist of much smaller crystalline subunits (10-. For instance, small- 
angle light scattering from CdS particles revealed that they were built of 
~ 50 A uniform globules, while electron microscopy studies substantiated 
the composite nature of uniform particles such as Ce02 (13) or Sn02 
(14,15). Finally, it has been shown with uniform colloidal copper oxide 
that the crystallite size in the final particles, determined by X-ray diffrac- 
tion, was the same as the diameter of the precursor subunits (singlets) 
(16). The latter finding added credence to an aggregation mechanism, 
rather than to a dissolution/reprecipitation process. 

The substructures have been identified also in particles of shapes 
different from spherical. For example, it was demonstrated that differ- 
ent morphologies of colloidal hematite (cK-Fe203) (12) were related to the 
natures of the precursor singlets (18-20). All this experimental evidence 
has led to the conclusion that the mechanism of the formation of uni- 
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form colloids is more complex than originally expected, and that it likely 
involves more than just a nucleation/diffusional growth process. 

Thus, the formation of monodispersed colloids poses two major the- 
oretical challenges. First, the morphology and shape selection of parti- 
cles formed by interplay of nucleation and aggregation processes must be 
explained. Second, the size- selection mechanism, i.e., the kinetics of gen- 
eration of narrow particle size distribution, must be identified. Several 
theoretical approaches utilizing thermodynamic and dynamical growth 
mechanisms (4,5,18,21-32) have been described in the literature. It has 
been established that models of aggregation of subunits can be devel- 
oped that yield a peaked and even sharpening with time particle size 
distribution (21-34). However, none of these attempts could fit quanti- 
tatively a broad range of experimental findings. In this work the first 
results of a new approach that explains the size selection are reported. 
Specifically, we argue that the growth by aggregation of subunits must 
be considered coupled to the process of formation of these subunits, via 
the time-dependence of the latter. 

Theories of coagulation of nanosize particles involve modeling of the 
particle transport towards each other and particle-particle interactions. 
The actual attachment event and possible detachment and restructuring 
processes, which presumably make the density of the final solids close to 



-5- 



that of the bulk material, must be then described. There exists exten- 
sive recent literature that generally addresses these and related issues in 
various situations. Let us first mention models of growth at surfaces (35- 
37) that typically assumed diffusional transport, or convective-diffusional 
transport, the latter being appropriate for flowing suspensions (38). The 
main ingredient of the simplest transport models is that the effects of 
particle-particle interactions are typically incorporated via average, phe- 
nomenological Boltzmann weights. Thus, the resulting deposit morphol- 
ogy, which depends on the precise way of particles fitting in and rear- 
ranging within the exiting structure, cannot be considered. Instead, the 
emphasis has been on the fluctuation properties of the topmost, outer 
growth layers and their scaling description (35-37). 

Traditional models of coagulation and nucleation usually assume dif- 
fusional transport and consider growth of particles and aggregates ei- 
ther via microscopic nucleation processes (21,39) or by particle-particle, 
particle-aggregate, and aggregate-aggregate adhesion on encounters (40,41) 
Models of sparse, percolative aggregates, typically predict the size distri- 
bution to peak at small sizes (42). More compact aggregate growth, such 
as in classical nucleation or coagulation, usually yields size distribution 
that grows with time, t, according to t^. Here z is termed the dynamical 
exponent and its values range from 0.2 to 0.5, depending on the model. 
Prescriptions for sharpening the size distribution by actually physically 
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removing particles of certain sizes in the course of the process, in com- 
peting nucleation and aggregation, were reviewed in (28). 

Recently, the simplest rearrangement processes after attachment, fol- 
lowing collision, have been considered, mainly by numerical modeling, 
including rolling to preferential sites and diffusional relaxation (43-48). 
However, rearrangement effects on the formation of monosize particles 
have not been explored. The experimentally observed role of aging in 
size-selection in fine particle synthesis will be addressed in the follow- 
ing sections. This role is not related to the restructuring of the growing 
particles. 

The rest of this article is organized as follows. In Section 2, a ki- 
netic model of secondary particle formation by singlet-capture dominated 
growth is considered. In Section 3, primary particle (singlet) formation 
by burst nucleation as well as other effects are incorporated to complete 
the model. Section 4 describes the gold sol used as the experimental sys- 
tem to test the model. Finally, Section 5 reports numerical simulations 
and comparisons with the experiments. 
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2. Kinetic Model of Particle Formation 

Uniform particle formation in colloid synthesis usually proceeds in 
several stages. In the initial induction stage solutes are formed to yield 
a supersaturated solution, leading to nucleation. The nuclei then grow 
by diffusive mechanism to form the primary particles (singlets), which in 
turn aggregate to form secondary particles. This latter process is some- 
times facilitated by changes in the chemical conditions in the system. For 
example, the ionic strength may increase, or the pH may change, causing 
the surface potential to approach the isoelectric point. Both effects result 
in a reduction of electrostatic barriers, thus promoting particle aggrega- 
tion. Formation of the final (secondary) particles, which can be of narrow 
size distribution, is clearly a diffusion-controlled process (2-7). Some of 
the various stages of the process can be more appropriately viewed as 
particle evolution rather system's dynamics. 

An important assumption of the present model is the experimentally 
observed property that the secondary particles are sufficiently sparsely 
positioned in solution to conclude that their evolution is largely inde- 
pendent. Furthermore, it is assumed that the particles are spherical, 
with density close to that of the bulk material, which is experimentally 
a very common case, and the focus will be on their size distribution due 
to growth by aggregation. These approximations will be discussed at the 
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end of this section. Several other simpHfying assumptions are made in or- 
der to zero-in on the essential ingredients of the theoretical modeling and 
make numerical simulations tractable. For instance, it is possible that 
singlets forming secondary particles may undergo restructuring and re- 
arrangements resulting in more compact structures and leading to shape 
selection. Such processes are neglected in the present model, but may be 
incorporated in future studies. 

Thus, let us focus on the two main stages in the process: the pro- 
duction of primary and secondary particles. The latter will be treated 
presently, while the formation of nuclei and their growth (aging) to pri- 
mary particles will be considered in the next section. Kinetic rate- 
equation-type models of aggregation typically utilize a master equation 
for the distribution of growing particles by their size. In the present 
model, the size will be defined by how many primary particles (singlets) 

were aggregated into each secondary particle, denoted by s = 1,2, 

The growing particles can adsorb or emit singlets and multiplets. The 
master equation is then quite standard to set up, as in, e.g., (33). In 
the present case, the process is experimentally documented to be highly 
irreversible, so that detachment can be disregarded. 

Furthermore, it is assumed that the singlets have diffusion constants 
larger than aggregates so that their attachment dominates the process; 
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this approximation will be discussed later. The master equation then 
takes the form 



= Ws-iNs_i - WsNs , for s > 1 . [1] 

Here Ns{t) is the time-dependent number density (per unit volume) of 
the secondary particles consisting of s primary particles. The attachment 
rate Ws will be taken from the Smoluchowski expression (49): 



Wsit) = 4TTRsDNi{t) , [2] 
where Rs is the radius of the s-size particle, given by 



= 1.2 rs^/^ . [3] 

The parameters r, the primary particle radius, and D, their diffusion 
constant, are experimentally available. The constant 1.2 in Eq. [3] was 
calculated as 



(0.58)-^/^ ~ 1.2 , [4] 

where 0.58 is the typical filling factor of the random loose packing of 
spheres (50). The fact that all primary particles were assumed to have 



the same radius r will be addressed in Section 3. Note that the rate in 
Eq. [2] involves the product RgD oc rD, where the latter, according to 
the Einstein formula for the diffusion constant, is not sensitive to the 
distribution of values of the radii r. 

In standard approaches to aggregation, the evolution of the popu- 
lation of singlets, which is not covered by Eq. [1], is obtained from the 
conservation of matter: 



which assumes that initially, at time t = 0, there are only singlets. In 
this approach, however, singlet-dominance cannot be assumed; specifi- 
cally, Eq. [1] and/or Eq. [5] need to be modified. Otherwise, the particle 
distribution is confined to small sizes, as confirmed by numerical simula- 
tions and other studies (46). Indeed, most singlets will simply combine 
into dimers, fewer trimers, etc., then the growth stops. The conven- 
tional approach has been to consider more general models, with discrete 
or continuous population balance, allowing multiplet-multiplet aggrega- 
tion, which adds terms in Eq. [1]; see, e.g., (21-34). These models will be 
commented on at the end of this section. Frequently, though, multiplet- 
multiplet aggregation models yield wide particle-size distributions, of the 
scaling type (33,35-37), and they cannot explain the size selection. 
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In this work we develop a different approach based on the observa- 
tion that the supply of "singlets," i.e., of primary particles, is in itself a 
dynamic process. Numerical calculations indicate that if the concentra- 
tion of the singlets were constant, i.e., if they were continuously generated 
to compensate for their depletion due to aggregation, then the resulting 
particle size distribution would be wide and peaked at small sizes, such 
that Ns ^ Nis~^^^ up to some growing cut-off s-value. However, if the 
availability of singlets is controlled to be a decaying function of time, then 
size selection can be obtained in some cases. For instance, when the rate 
p{t) at which the primary particles are formed (per unit volume) was cho- 
sen to be decaying according to a power-law, then preliminary numerical 
calculations yielded a single-hump size distribution for the secondary par- 
ticles. Note that the equation for Ni{t) must be modified by replacing 
Eq. [5] with the relation 



with the initial values -/Vs(O) = for all s = 1, 2, 3, 

Equations of the type considered here, with singlet-capture domi- 
nance of the dynamics and several ad hoc singlet-input rate functions, 
have been described in the literature (51). The emphasis of these studies 
has been on cases which are exactly solvable, for instance, one-dimensional 
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versions, and those where the size distribution shows self-similar "scaling" 
behavior. The goal of the present work is quite different, i.e., to actually 
model the primary-particle input rate p{t) by using the burst-nucleation 
approach; see the next section. 

Let us now address some of the approximations involved in our treat- 
ment of the aggregation process. These include, for instance, ignoring 
multiplet mobility and multiplet-miltiplet collisions, as well as the effect 
of mobility of aggregates (multiplets) on the diffusion constants used in 
the rate expressions, etc., especially in the beginning of the process when 
most aggregates are small. There are established methods in the liter- 
ature that avoid many of such difficulties (21-34) and, in fact, some of 
the models for various reaction rates lead to particle-size distributions 
peaked and even sharpening with time. However, the main point of our 
work, which we believe is new and crucial to obtaining narrow size dis- 
tributions, has been that the growth of the secondary particles must by 
coupled, via the rate of generation of singlets (primary particles), to the 
time-dependence of the process of formation of the latter; see the follow- 
ing sections. 

Thus, we intentionally took the simplest possible models of both pro- 
cesses, the primary and secondary particle formation. This choice has the 
following advantages: it simplifies numerical simulations and thus allows 
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to scan a wider range of model parameters; it avoids introduction of un- 
known microscopic parameters. As a result, for instance, our description 
of the aggregation process in this section has no adjustable parameters; 
they are all experimentally available. 

Similarly in the following sections, for primary particle formation, 
we only utilize one adjustable parameter, the effective surface tension, 
and even that parameter turns out to be close to the experimental bulk 
value. However, we recognize that more sophisticated modeling can im- 
prove consistency with experiment, perhaps at the expense of additional 
assumptions and parameters, and we intend to explore this avenue of 
investigation in future studies. In fact, we carried out preliminary nu- 
merical simulations allowing for the dimer diffusion and attachment to 
larger aggregates and showing trend of improved consistency with the 
experiment. 

We furthermore note that the approximation of restricting the aggre- 
gation process to only the smallest particles sticking to the larger particles 
has been widely used in the literature, for instance, (22,27). Considera- 
tions of colloid stability have been typically used to justify such approxi- 
mations, and we note that detailed arguments of this sort would require 
additional microscopic parameters in the model. 

We also comment that generally in diffusion-limited growth the ag- 
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gregates are expected to be fractal (35-37,52). In this work we avoid the 
issue of shape selection and internal structure of the secondary particles; 
see a review (21). There is obviously some restructuring going on dur- 
ing the secondary particle growth for the experimental system described 
in Section 4, because the particles are clearly spherical throughout the 
growth process. Actually, in some other experiments, e.g., (14), the sur- 
face of the secondary particles is initially "hairy" and it gets smoother 
at later times. In our case the internal restructuring processes are fast 
enough so that the shape/morphology of the aggregates is always spher- 
ical/compact. One could propose that the effects of restructuring in our 
model make the coefficient in Eq. [3] an adjustable parameter rather than 
a fixed number. We did not explore this matter. 
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3. Modeling the Primary Particle Formation 

Primary-particle production rate p{t) will be calculated by using 
the relations available in the literature (8,9,34,53) for the fast nucleation 
process. Such thermodynamic models of "burst" nucleation of more or 
less uniform dispersions have been originally formulated in (8,9). Model- 
ing the rate of formation of primary particles (singlets) requires setting 
up a master equation, where the rate of growth is determined by the 
Boltzmann factor with the thermodynamic free energy difference AG, 
multiplied by —l/{kT), in the exponent. This approach in turn requires 
modeling of the free energy of the growing embryos (sub-critical nuclei); 
in the simplest model one can use generic volume-plus-surface energy 
forms. 

Let us term "solutes" the species (atoms, ions, molecules) which 
serve as monomers for the primary-particle nucleation. Given the con- 
centration c{t) of solutes, which is larger than the equilibrium "satura- 
tion" concentration cq and approaches cq for large times t, the rate of 
formation of critical nuclei can be written as (34,53) 

p(t) = Wj/^Oc^e-^^-/*^ , [7] 

which is based on the diffusional capture of solutes, whose effective radius 
is denoted by a, diffusion constant by T>, and n is the number of solutes 
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in an embryo. The subscript ens refers to values calculated at the critical 
nucleus size. Note that c{t) = ni{t). 

The expression [7] involves the following ingredients. For embryos 
of size n < ricns, one assumes that solutes can be captured and emit- 
ted fast enough so that the size distribution is given by the equilibrium 
form. Thus, the factor c^-^Gcns/kT -^^ j-yj fQ^^Q^g from the assump- 
tion that sizes up to ricns are "thermodynamically" distributed, according 
the Boltzmann form. For sizes larger than ricns the dynamics is usually 
assumed to be fully irreversible, corresponding to an unbound growth 

1 /3 

by the capture of solutes. The factor A-KaricnsT^c in Eq. [7] is thus the 
appropriate version of the Smoluchowski growth rate similar to that in 
Eqs. [2,3]. The filling-fraction correction factor was absorbed in the defi- 
nition of a to simplify the notation; it will be specified in Section 5. 

For the free energy of the n-solute embryo, the known expression is 
used: 

AG = -nkT In (c/cq) + A-KO^n^f^a , [8] 

which involves the bulk term, proportional to n, and the surface term. 
The standard form of the bulk term was derived as follows. It is assumed 
that the entropic part of the free-energy change between the solid and 
solution phases can be calculated as the entropy in the supersaturated 



liquid suspension of solutes of concentration c, as given by the dilute 
(noninteracting) expression of the "entropy of mixing," defined, e.g., in 
(54). The surface term in Eq. [8] corresponds to the assumption that 
the growing embryos are spherical, of radius an^/^, and introduces their 
effective surface tension cr, which is usually assumed to be comparable to 
the bulk surface tension. It will be shown later that in the present case 
the results are very sensitive to the value of cr. Obviously, all the above 
expressions only apply for large n. It has been a common practice in the 
literature to use them for all n, as one of the approximations involved in 
a model. 

It follows that both 

TLcns cind ^Gcns ^^rc expHcit functions of c(t). 



SkTln (c/co). 



[9] 



2567r3aV3 
27(/cr)2[ln(c/co)]' 



^Gcns = ,2 ' [10] 



where the critical value ricns was calculated from dG/dn = 0. 

As the next step, the decrease in the concentration of solutes owing to 
the formation of critical nuclei is taken into account. Usually, one accepts 
that for n > ricns the primary particles grow (age) largely by absorbing 
diffusing solutes, and as in the preceding section we ignore here more 



complicated processes such as capture of small embryos, dissolution, etc. 
Simultaneously, the primary particles are also captured by the secondary 
particles. In the present model, it is assumed for simplicity that the 
primary particles are captured fast enough by the growing secondary 
particles so that the effect of their aging on the concentration of solutes 
can be ignored. Furthermore, it has been generally recognized that aging, 
when significant, tends to sharpen the size distribution (21,34). Thus, 
the primary particle radius r, introduced in the preceding section, will be 
assumed to have a single, experimentally determined value, although in 
reality they have a finite-width, albeit narrow, size distribution. These 
are among several simplifications on which the present model is based; 
they can be relaxed in later studies. 

Accordingly, one can write 



^ = -ncnsp{t) , [11] 

which means that the concentration of solutes is "lost" solely due to the 
irreversible formation of the critical-size nuclei. Collecting all the above 
expressions, one gets the following equations for c{t) and p{t): 



dc _ 1638471^ a^a'^Vc^ \ 2567T^a^a^ . 

di~~81{kT)^ln{c/co)f'''''^[27{kT)^ln{c/co)f ^ ' 
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Pit) 



SkTln (c/co) 



exp 



27(/cr)3[ln(c/co)]' 



[13] 



It should be noted that replacing the distribution of the primary 
particle sizes by a single, experimentally measured average value of r, 
violates the conservation of matter. Thus, in the present model only the 
shape of the secondary particle size distribution is relevant. The absolute 
number densities Ng must be rescaled to correspond to the actual amount 
of the solid matter per unit volume. The latter data are usually available 
experiment ally. 
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4. Experimental System 

In order to test the model developed in Sections 2 and 3, disper- 
sions of spherical gold particles were used, which were produced by the 
reduction of chloroauric acid (HAuCU) with ascorbic acid (55). The sim- 
plicity of the chemical reactions involved, the stability of the solid phase 
formed, and the possibility to either measure or estimate all the necessary 
parameters make this system suitable for testing the described model. 

Several methods of preparation of stable monodispersed gold sols by 
the reduction process have been reported (56-58). Most of these studies 
have been focused on approaches that avoided the aggregation process 
(i.e., very dilute and/or heavily stabilized systems). In contrast, in con- 
centrated dispersions, the aggregation process is assured, resulting in the 
formation of particles consisting of a large number of subunits. The pre- 
cipitation procedure used in the described study (55) has resulted in large 
spherical gold particles of a narrow size distribution. After a short induc- 
tion period (up to 6-8 sec) the nucleation occurs, followed by immediate 
aggregation. The total reaction time varies from 3 to 20 sec, depending on 
the experimental conditions selected. Scanning electron micrographs in 
Figure 1 illustrate the resulting particles, at two different magnifications, 
obtained by rapid addition of chloroauric acid to a solution of ascorbic 
acid under intense agitation. 
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The field emission microscopy, Figure 2, clearly reveals the presence 
of the subunits, having an approximate size of 30 to 40 nm. The cal- 
culated value for the packing fraction of the subunits in the aggregated 
particles, ~ 58%, was based (50) on the experimental tapped density of 
dry powder. This value is characteristic of a random loose packing, usu- 
ally expected in the formation of rapidly assembled systems (50). The size 
of the primary particles can be also estimated from the X-ray diffraction 
measurements by the Scherrer formula, if one assumes that the subunits 
are perfect crystals. Several calculations, using different peaks, have gen- 
erated values between 30 and 42 nm, which is in excellent agreement with 
the data from electron microscopy. Finally, in an attempt to prepare and 
characterize the primary gold particles before undergoing aggregation, the 
same precipitation process was carried out in very dilute solutions where 
the metal particles were electrostatically stabilized. Figure 3 shows the 
size distribution of these particles, as determined by dynamic light scat- 
tering. Again, the modal diameter of ~ 40 nm agrees with the previous 
values for the precursor subunits. 
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5. Results and Discussion 

Numerical simulations required to follow the time evolution of the 
kinetic equations turned out to be large-scale. Therefore, the present 
testing of the model has been restricted to one randomly selected set of 
experimental parameters. However, we actually varied numerically all the 
parameter values and found that the calculation results are affected to 
various degree by them. The parameter values will be discussed roughly in 
the order of increasing sensitivity of the numerical results. For consistency 
all input data are in the MKS system of units. 

The radius of the primary particles, 

r = 2.10- 10~^m , [14] 

was obtained as described in the preceding section, which is within the 
range of 0.5 • 10^^ m to 5 • 10~^ m, typical for the system under consid- 
eration. The value in Eq. [14] applies to the experiment for which the 
initial concentration c(0) was 

c(0) =6.0-10^^m-^ ; [15] 

c(0) was calculated from the concentration of the gold solution used in 
the preparation of the dispersion, and it yields the initial condition in 
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Eqs. [12,13]. The diffusion constant D of the primary particles was ob- 
tained from the Einstein formula, 



D = 1.03- IQ-^V^sec"^ , [16] 

with 

kT = 4.04 • 10-2^ J . [17] 

The saturation concentration of gold in solution, cq, is not well known 
and is expected to depend somewhat on the experimental conditions. 
Using 2 • 10~^^ moldm"^ (59) yields 

Co ~ 1 • 10^^ m-^ . [18] 

The results for the particle size distribution are not particularly sensitive 
to this parameter, because it enters under logarithm in Eqs. [12,13]. The 
solute diffusion constant, 

r> = 1.5 • IQ-^m^sec-^ , [19] 

was estimated similarly to D in Eq. [16], using the Einstein formula, with 
the radius of the gold atom of 1.44 • 10~^°m (60,61). The applicability 



of this formula to single atoms may not be exact, but the particle size 
distribution is not too sensitive to this parameter value: a decrease in V 
shifts the calculated distribution to somewhat smaller aggregate sizes. 

It was established numerically that the size distribution of the sec- 
ondary particles was sensitive to the values of a, the effective atomic 
(solute) radius, and to the surface tension a. Note that a was defined to 
relate the number of solutes n in a growing primary particle to its radius, 
given by am}^^ . It is assumed that the primary particles are largely crys- 
talline; thus, the best choice of a is such that 47ra'^/3 is the volume per 
atom, including the attributable part of the surrounding void volume, in 
bulk gold. Consequently, 

a= 1.59- 10-^° m [20] 

was obtained by dividing the radius of the gold atom (1.44 • 10~^°m) 
by the cubic root of the volume filling fraction, 0.74, of the crystalline 
structure of gold (50). 

The effective surface tension of nanosize gold embryos in solution, cr, 
profoundly affects the numerical results. Unfortunately, even the bulk- 
gold value, which is of order 



(J ~ 0.58 to 1.02 Nm 



-1 



[21] 



is not well known (62), and it may differ from that of the "nanosize" 
solids. Given this fact, a was chosen as the only adjustable parameter in 
the model. 

Experimentally, the time scale on which the secondary particle growth 
effectively terminated was about 8 to 10 sec, which does not include the 
"induction" stage. In Figures 4 though 6, the results of the numeri- 
cal simulations of the kinetic equations are presented with parameters 
as specified above, for three different values of cr, which clearly demon- 
strate the sensitivity to the choice of this parameter. In Figure 4, the 
case (7 = 0.51 N/m illustrates growth that already reached saturation for 
times up to 10 sec. It should be noted that the distribution evolves quite 
slowly with time. Initially, it is heavily weighed in the small-aggregate 
regime. Later on, the large-size peak develops and eventually dominates 
the distribution. By varying a near the expected range, given in Eq. [21], 
it was found that, for times up to 10 sec, all a values yielded smaller 
average sizes than the experimentally measured one. 



Seeking a that would give the largest secondary-particle size resulted in 




[22] 



Wfiued = 0-57 ± 0.04 Nm-1, 



[23] 



which agrees well with the bulk value in Eq. [21]. 

Figure 5 shows the size distribution for a = 0.57 N/m. The growth 
did not reach full saturation on the relevant time scales, and the peak par- 
ticle radius at t = 10 sec, of Rg — 0.32 /im, is somewhat smaller than the 
experimental value in Eq. [22]. The width of the distribution, of ~10 %, is 
close to that established experimentally. Considering the approximations 
involved in the model, only semiquantitative agreement with the exper- 
imental data should be expected. Since the key feature of the model 
is the prediction of the narrow-width distribution of secondary particle 
sizes, the overall consistency with the experimental results is gratifying. 

As the value of a is increased, the large-size peak does not fully de- 
velop on the relevant time scales, as exemplified by the case a = 0.63 N/m 
shown in Figure 6. The reader should be reminded that, owing to the 
absence of the conservation of matter in this model, the number densities 
Ng must be rescaled according to the actual amount of the solid matter 
per unit volume, if the comparison of the calculated and experimental 
distribution is attempted. This refinement will be addressed in future 
studies. 

Asymptotically, the particle-size distribution "freezes" for large times. 



i.e., the particle growth actually stops in this model as opposed to the 
"scaling" growth studied, for instance, in (51). To demonstrate this prop- 
erty, note that Eq. [11], with Eq. [9], can be integrated in closed form to 
yield 

/*'.-=(i^)"-K?)-(f)]' 



where 



F{x) = x (In xf -S{\nxf -\-6{\nx) -Q . [25] 

The left-hand side of Eq. [24] is the total number of primary particles 
produced by the time t. Obviously, from Eqs. [24,25] this number is 
finite as t ^ oo, i.e., in the limit c{t) cq. The supply of the primary 
particles, manifested by the peak at small sizes for short times, in Figures 
4-6, is essential for the large-size peak in the distribution to develop and 
grow, because the present model assumes the growth of the secondary 
particles to be solely by singlet capture. 

An interesting mathematical construction is suggested by ideas de- 
veloped in (51). Let us consider the quantity r defined by 
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t 

r = J Ni{t')dt' . [26] 



If the independent variable is changed from t to r in Eq. [1-4], the result- 
ing equations are linear in Ns>i. One can then easily show that the only 
way to have a normalizable stationary size-distribution as t ^ oo is to 
have r(t) approach a finite value for large t. This quantity was calculated 
numerically for the a values used in Figures 4-6, and the results, shown 
in Figure 7, confirm the earlier observations. Specifically, the growth pro- 
cess saturates fast for a = 0.51 N/m. For the two larger a values there is 
still some variation for the time scales of order 1 to 10 sec. This function 
is useful in identifying the time scales of the growth process. 

In summary, a model is proposed for the synthesis of particles with 
the size distribution that is peaked at an average value corresponding to 
a large number of primary particles in a final secondary particle. For the 
experimental gold-sol system the model has worked reasonably well: the 
average size, the width of the distribution, the time scale of the process, 
and even the fitted effective surface tension were all semiquantitatively 
consistent with the measured, known, or expected values. The present 
model has involved numerous simplifying assumptions. Further studies 
are needed to incorporate additional effects in the model and test it for 
a wide range of experimental systems. The main conclusion has been 
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that multistage growth models can yield size-selection as a kinetic phe- 
nomenon, which has been observed in a large number of experimental 
systems. 
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Figure Captions 

Figure 1: Scanning electron micrographs, at two different magnifications, 
of final (secondary) gold particles obtained by adding rapidly 100 cm^ 
of an aqueous solution of HAuCU (0.5 mol dm~"^) into 400 cm"^ of an 
aqueous solution of ascorbic acid (0.5moldm~^), at room tempera- 
ture; see (50). 

Figure 2: (A) Field emission microscopy image of gold particles shown 
in Figure 1. (B) Enlarged image of the darkened area in (A). 

Figure 3: Particle size distribution of primary particles (singlets) of gold 
prepared in a dilute system (2 • 10~^ moldm"^ HAuCU), determined 
by dynamic light scattering. 

Figure 4: Distribution of the secondary particles by their sizes, cal- 
culated using cr = 0.51 N/m, for times t = 0.001, 0.01, 0.1, 1, 10 
sec. Parameter values used in the numerical calculation are given in 
Section 5. 

Figure 5: The same plot as in Figure 4, using a = 0.57 N/m, for times 
t = 0.1, 1, 10 sec. 

Figure 6: The same plot as in Figure 4, using a = 0.63 N/m, for times 
t = 0.1, 1, 10 sec. 

Figure 7: The function r(t), defined in Section 5, calculated for the 
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values of a corresponding to Figures 4-6. 
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